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The neutrino community is convinced that if the LMA solution is indeed correct, the prospects 
for next-generation, long-baseline neutrino experiments are very exciting. In this work we will argue 
that the prospects are not less exiting if the answer to the solar neutrino puzzle lies in the LOW 
region. For this purpose we explore the oscillation probabilities which are accessible to experiments 
which employ conventional neutrino beams. We consider the electron (anti)neutrino appearance 
channel, which provides information regarding not only the small C/e3-element of the leptonic mixing 
matrix but also the neutrino mass hierarchy, and also include information, which should come from 
. . . the muon (anti)neutrino disappearance channel, regarding the magnitude of the atmospheric mass- 

I squared difference (|Ami3|) and mixing angle {O23). We comment on the presence of ambiguities 

■ and the experiments needed to solve them. 
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p ,i I. INTRODUCTION 

(D ; 

' There is unambiguous evidence of new physics in the leptonic sector of the standard model. Furthermore, all 
, experimental neutrino data ^, ^, |^] are consistent with the hypothesis that neutrinos have mass and mix. This 
• being the case, the goal of this and the next generation of neutrino experiments is to determine the neutrino masses 
[ and the leptonic mixing matrix to the best of their capabilities. 
^ . One of the issues which is to be resolved by the present generation of neutrino experiments involves the determination 
' of the so-called solar parameters, Amf2 s-nd tan^ di2. The current data strongly require tan^ 612 to be of order unity, 
^ ; while 10"^ eV^ < Amfj ^ 10"^ eV^ (the large mixing angle solution (LMA)) or Aml^ < lO""^ eV^ |, |, 0. The 
I ' latter will be referred to in this paper as the LOW solution.* While global analyses indicate that the LMA solution 
0^ , is significantly preferred, it is fair to say that the LOW solution still remains a realistic possibility 

■ In any case, the KamLAND reactor neutrino experiment ^ will, in the near future, definitively establish whether 
the LMA solution is driving solar neutrino oscillations. If KamLAND observes an oscillation signal, it will not only 
exclusively select the LMA solution but also "pin-point" the value of solar parameters |ll| . If KamLAND does 

' not observe an oscillation signal, the LMA solution will be unambiguously (modulo CPT violation |l^) ruled out pl , 
and we will have to wait for the Borexino solar neutrino experiment |13| , which is scheduled to start data-taking in 
2003, in order to definitively piece the solar neutrino puzzle p^ . 
D , If the LMA solution is indeed correct, the prospects for next-generation, long-baseline neutrino experiments are 
1-^ ' very exciting. Not only can these experiments probe the remaining mixing angle of the leptonic mixing matrix, but 
J> ] there remains the possibility of also observing CP-violating effects. This scenario has been thoroughly studied, for 

' different long-baseline experiments (superbeams, neutrino factories, etc), in the literature |l7| . 
/\ ' In this work we will explore the physics capabilities of future long-baseline oscillation experiments if the LMA 
^ , solution is not correct. We will argue that a rich physics program lies ahead despite the fact that the door to leptonic 
■ - - ' CP violation will be closed (indeed, the study of any "solar effect" in terrestrial oscillation experiments becomes 
virtually impossible). First, we discuss the oscillation probabilities which are accessible to experiments which employ 
conventional neutrino beams (which are mostly composed of muon-type (anti)neutrinos). We consider both the 
electron (anti)ncutrino appearance channel, which provides information regarding not only the small [/e3-element of 
the leptonic mixing matrix but also the neutrino mass hierarchy, and the muon (anti)neutrino disappearance channel, 
which provides information regarding the magnitude of the atmospheric mass-squared difference (jAm^gl) and mixing 
angle (023)- A comparison of both expressions reveals that the information provided by the disappearance channel 
leads to "ambiguities" ||l^ in the determination of the currently unknown oscillation parameters. Far from being a 
nuisance, these degeneracies are an indication that more physical quantities can be measured given the correct amount 
of information. 

Next, we simulate and analyse "data" from a finely segmented iron-scintillator detector located off-axis from the 
future NuMI beam at Fermilab. In measuring the oscillation parameters, we include the realistic uncertainties on all 
relevant parameters, and consider several scenarios in order to fully study all degeneracy issues. We find that one 
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* The current data allow, at th<l.^''tbrcc sigma" level, not only the LMA solution, but also the LOW, QVO, and VAC solutions, depending 
on the data analysis. Sec, M, fel [7| for details. 



2 



should be able to determine the neutrino mass-hierarchy by combining information from neutrino and antineutrino 
oscillations if |J7e3p is large enough (solving one "degeneracy issue"), but not whether ^23 is less than or greater 
than 7r/4. We comment on how this degeneracy might be attacked. Finally, we analyse the issue of how well the 
atmospheric mass-squared difference should be known before choosing the location of the off-axis detector. 



II. THE DISAPPEARANCE CHANNEL 



We are interested in baselines around 10'^ km and order 1 GeV neutrinos energies, chosen in order to explore the 
first oscillation peak. We hope to show, contrary to some claims, that there is no need to go beyond the first peak in 
order to disentangle the oscillation features. For this class of experiments, the parameter a (see Appendix A), which 
measures the strength of matter effects, is given by 



l^flGpn^E 2.8 10-4 / E 



Amf^ Aml^ V GeV 



(II.l) 



where rig is the electron number density in the Earth and E is the energy of the neutrino (see Appendix A for 
definitions and conventions). 

In the region of interest, a is relatively small (a ~ .1 — .2) such that, up to quartic terms in |J7e3| = Sia, the muon 
neutrino survival probability can be written as (complete formulas are given in Appendix A) 

Pi^.-^i^^.) = 1 - 4^2344 sin2[Ai3]-s23^^1f^sin2[Ai3(l- a)]. (II.2) 

Note that we have also expanded A13 — Aisjl — a|(l + 0{a\Ue3\'^))- Matter effects, therefore, only modify the 
vacuum survival probability by terms which are 0(|[/e3pa), and can safely be neglected. The vacuum expression can 
be written as 

P{iy^ ^ J^m) - ^('^M ^m) = 1 - 4|C/^3p(l - |C/M3p) sin2(Ai3); (II.3) 
I [/^3 1 = sin 6*23 cos 6*13, (II.4) 

meaning that one can measure jATO^gj (from the oscillation frequency) and sin^ 2^23 = 4|C/p3p(l — |?7^3p) (from the 
magnitude of the oscillation) . For | Ues \ <^ 1 



sin^ 29^^ = sin^ 26I23 



sm 26*23 



(11.5) 



In the limit of vanishing \Ue3\, sin^ 26*23 — sin^2023, while a precise enough extraction of sin^ 26*23 from sin^ 2(?23 
is sensitive to the value of |l7e3|- Indeed, the fact that |C/e3p is not known induces an irreducible uncertainty on 
extracting sin^ 2^23 from z/^ disappearance experiments as large as A(sin^ 2^23) ^ 0.02, depending on the value of 
sin^ 2023 and the uncertainty on |C/e3p-^ 

The information needed in order to interpret the appearance channel (as will be discussed in the next section) 
is sin^023 (not sin^ 26*23) and (Am^g (not |A 

'^131)' This implies that, as far as the appearance channel is concerned, 
the information from the disappearance channel is four-fold degenerate, meaning that a unique measurement of 
(sin^26l^f, |Am^3|) "splits" into 

(sin^ 923, Ami,) ^1 - sin^ 26^^f ) , ±|Am?3|^ , (II.6) 

where the subleading |?7e3p effects have been ignored. 

This parameter degeneracy is depicted in Fig. |l|, where we present the one and three sigma confidence level (CL) 
contours in the (sin^ ^23 x |C^e3p) and (sin^ 6*23 x ATO^3)-planes, assuming that one has measured |ATOf3| = (3±0.1) x 
10-3 eV^ and sin^ 26*^1 = 0.91 ± 0.01 (to be exact, the input point corresponds to sin^ 6*23 = 0.35 and |f7e3p = 0.01, 
as indicated by the (blue) stars in Fig. |l|). We also impose the CHOOZ bound jl^ as a hard cutoff, |t/e3p < 0.05, 



t It is curious to note that for "maximal mixing" in the atmospheric sector (623 = ""Z^), there are only |C/c3|^ corrections. 
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and that the measurements of the effective atmospheric angle and the absolute value of the atmospheric mass-squared 
difference are uncorrelated. This is a simplifying assumption, which would be a good approximation if future long- 
baseline disappearance experiments run at distances and energies which are around the minimum of the survival 
probability. In this case, and assuming a good energy resolution and calibration, the mass-squared difference can be 
obtained by the position of the minimum, while the sin^ 2^§f is obtained by the magnitude of the suppression. The 
capabilities of long-baseline experiments to measure the atmospheric parameters is currently under investigation , 
and clearly depend on the detector-type and the neutrino beam. 
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FIG. 1: One and three sigma confidence level contours in the (sin^ 623 x |[/e3p) and (sin^ 623 x Ami3)-planes, obtained from 
the "measurement" lAmfal = (3 ± 0.1) x 10~^ eV^ and sin^ 2623 = 0.91 ± 0.01, depicted in the inset located in the center. See 
text for details. 

In each plot of Fig. [l], the CLs are defined for two degrees of freedom, i.e., the absent parameter in each figure has 
been "integrated out." On the bottom plot, we present the confidence ellipses for different fixed values of |C/e3p- The 
appropriate one and three sigma level regions are to be interpreted as the envelope of the different inner and outer 
ellipses, respectively. 

Several comments are in order. First, the alluded to four- fold degeneracy is clearly visible in the bottom plot of 
Fig. 0. Second, the fact that |J7e3p is not known implies a significant increase of A(sin^023), as can be observed by 
comparing each individual ellipse in the bottom plot with the size of the envelope. The situation will improve once 
more information (from the appearance channel) regarding the magnitude of Ue3 is included. 

To conclude this section, we point out the simple but important fact that a Gaussian measurement of sin^ 2^23 
does not translate into a Gaussian measurement of sin^^23. This is true even when |J7e3| = 0. This can be easily 
understood in the following way. If the measurement of sin^ 29 is Gaussian, it can be transformed into a function 
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(sin^ 26 — xY jior^Y 1 where x = sin^ 20* is the measured value of sin^ Q with error ^Ox- This translates into 

where yi,2 are the solutions to y(l — y) — a;/4 = 0. If yi,2 are separated enough (in units of cr^) the splits into the 
sum of two parabola, with error Uy — (Jx/\^{vi — 2/2) | (of course, \yi — y2\ — \ cos 26**1). This will not be the case of 
maximal mixing {yi = 2/2 = 0.5), where the function will be a quartic monomial which cannot be approximated by 
a parabola. This implies that a measurement of sin^ 29—1 with an uncertainty of 0.01 will imply sin^ 6 = 0.5 ± 0.05, 
but the ±0.05 one sigma uncertainty is very non-Gaussian, meaning that, for example, the three sigma error bar is 
not ±0.15 (as a matter of fact, the three sigma error bar is significantly smaller than ±0.15). 

III. THE APPEARANCE CHANNEL 

The electron neutrino appearance probability is given by (again, for small a) 

P i^, ^ ^e) = 4377^^ sin^ [Ai3(l - «)] , (IIL8) 



while the approximate expression for antincutrinos is given by Eq. (III. 8) if one exchanges a ~a. (For exact 



expressions see Appendix A). Note that Eq. ( III. 8 ) can be significantly modified by the presence of matter (at 0{a)). 

The dominant effect will arise from the — modification of the overall magnitude of the oscillation probability, 
since we assume that the available experimental information will not be able to properly "see" the position of the 
oscillation maximum. This implies that the Ve appearance rate for a fixed value of |C/e3p is enhanced [suppressed] 
with respect to the pure vacuum one if the neutrino mass hierarchy is normal (a > 0) [inverted, {a < 0)]. The 
opposite happens in the antineutrino channel. If the mass hierarchy is not known, a degeneracy in the determination 
of |C/e3|^ surfaces, as a fixed value of |?7e3|^ and a normal mass hierarchy will yield the same number of appearance 
events within some energy window as a larger value of [J/esj^ and an inverted mass hierarchy. This ambiguity can 
be solved by, for example, combining the information obtained with the neutrino and antineutrino channels (whose 
matter effects are opposite) or by changing the baseline of the experiment, such that the significance of the matter 
effect is modified. Note that precise information regarding the location (i.e., the incoming neutrino energy) of the 
minimum of the survival probability and the maximum of the appearance probability will also identify the 
neutrino mass hierarchy. This information, however, is hard to obtain in the appearance channel. 



Another degeneracy ensues if sin^ 6*23 7^ 0.5. This comes from the fact that Eq. ( 111.8 ) is directly proportional to 
sin^ 023, while, as discussed in Sec. ^ (see Fig. |l|) the disappearance channel cannot distinguish sin^ 6*23 from 1 — sin^ ^23 
(ignoring small jJ/esj^ corrections). This means that for every (|f7e3|^, sin^ 6*23) pair there is another choice for the 
oscillation parameters (namely, (|?7e3|^ tan^ ^23, cos^ ^23)) whic h yield s exactly the same P(f^ — > Ve) as a function of 



energy, as long as the approximations that go into writing Eq. ( IIL8| ) apply. Contrary to the "hierarchy degeneracy," 
the degeneracy in the atmospheric angle cannot be solved by comparing conjugated channels, as it is not modified 
by the presence of matter. Similar experiments at different baselines and neutrino energies will also have no effect. 
In order to eliminate this degeneracy it is necessary to either look at different oscillation modes (which come with 
distinct ^23 dependences) or to probe P{vi_i — > i^e) around a = 1, where the value of |J7e3|^ affects the appearance 
rate in a less trivial way. 

Candidates for different oscillation modes include the search for the disappearance of electron-type neutrinos or 
i^e ^ oscillations. For |?7e3|^ ^ 1 

P(iye->i^e) = 1- — ^^sin2[A23(l-a)], (IIL9) 
(1 - a) 

P {v, ^ vr) - c^g-i^ sin^ [A23(l - a)] , (III.IO) 

(1 - ay 

which depend only on jJ/esj^ or on cos^ 023|C^e3|^- We will comment more on this issue in the next section. 

IV. SIMULATIONS OF THE OFF-AXIS EXPERIMENT(S) 

We simulate and analyse "data" for electron (anti)neutrino appearance in a highly segmented iron-scintillator 
detector located off-axis from the NuMI neutrino beam. The detector is described in |M], together with the off-axis 
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neutrino beams and the reconstruction efficiencies, and we refer readers to it for more details. We will always assume 
that the data consists of 120 kton- years of running with a "neutrino beam" (see [|l6|) plus, whenever applicable, 300 
kton- years of running with an "antineutrino beam." We assume that the detector is located 12 km off- axis and 900 km 
away from the neutrino source (see jl^), unless otherwise noted. We further include in the "data" analysis the fact 
that jAm^gl is measured with a precision of A(|Am^3|) — 0.1 x 10~^ eV^, while sin^ 262^ is measured with a precision 
A(sin2 26*^1) = 0.01. This type of precision has been quoted, for example, in studies of the physics capabilities of 
a future JHF to SuperKamiokande neutrino program and we assume that the study of muon disappearance at 
the off-axis detector being considered here should yield similar precision. Both errors are considered to be one sigma, 
uncorrelated, Gaussian errors. -I- 

In order to discuss all relevant physics issues, we will consider several scenarios: sin 623 — 0.35,0.5,0.65, Aml^ = 
±(2, 3, 4) X 10"'^ eV^ and |C/e3p < 0.05. These choices are meant to be close to the current best fit of the atmospheric 
neutrino data [tI or at the boundaries of the currently allowed 99% CL region. For concreteness, we have fixed 
Am^2 = 0.93 X 10 eV^, sin^ 812 = 0.33. It should be clear that any other LOW value for the solar parameters would 
have yielded identical results. 

First, we determine the capability of the setup in question to observe a signal. We do that by computing the 
sensitivity to appearance as a function of the relevant oscillation parameters. The three sigma sensitivity CL curve 
for 120 kton-years of neutrino-beam running (see ||l6| for details) is depicted in Fig. |2|in the (|J7e3p x Am^gj-plane, for 
both neutrino mass hierarchies (signs of Am^g) and three different values of sin^ 623- By sensitivity we mean that, if 
|?7e3p, Amjg correspond a point in the curve, the probability that the number of events observed after 120 kton-years 
of running is due to background is 0.27%. 
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FIG. 2: Three sigma confidence level sensitivity of 120 kton-years of "neutrino beam" data, 12 km ofT-axis and 90 km away 
from the NuMI beam, as a function of Ami3 and |(7e3p, for sin^ 623 = 0.35, 0.5, 0.65. The hatched region is currently ruled out 
by the atmospheric neutrino data at the three sigma confidence level. 



As discussed in Sec. Ill, the sensitivity is greater for larger values of sin 623 and for the normal neutrino mass 
hierarchy, due to a larger number of I'e-appearance candidates. Not surprisingly. Fig. ^ indicates that the optimal 
sensitivity occurs for Am^g ~ 3 x 10~^ eV^, as we have chosen the neutrino-beam energy profile {i.e., the off-axis 
distance) and experimental baseline in order to obtain A13 ~ 7r/2 for this particular value of the mass-squared 
difference. 



Once a signal is detected, the next step is to try and determine the neutrino mixing parameters. 



and sin 623. We do this for several different input values. The two-sigma confidence level measurements in the 



As noted in Sec. M, a Gaussian error in sin^ 2023 does not leads to a Gaussian error for sin^ 623- This is appropriately taken into account. 
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(Am^3 ^ l^^esp), (|f^e3p X sin^023), and (Am^o x sin^ 023)-planes obtained after 120 kton-years of "neutrino-beam 
data" are depicted in Figs. |3[ ^, and H, for sin 6 = 0.5, 0.35, and 0.65, respectively, Am^g = +3 x 10^"^ eV^ and 
l^^eSp = 0.008. In the case of "maximal" atmospheric mixing (Fig. H), there are still two disconnected solutions, 
centered around two distinct values of Am^g and |C/e3p- In this case, one would obtain at the two-sigma level 
0.005 < |C/e3p < 0.02. For nonmaximal atmospheric mixing (Figs. ^, the situation is even more confusing, and 
there are four disconnected solutions centered at distinct values of all the three parameters of interest. 
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FIG. 3: Two-sigma confidence level regions in the (Am^s x |C/e3|^), (|C/e3l^ ^ sin^623), and (Ami3 x sin^ 6'23)-planes obtained 
after 120 kton-years of "neutrino-beam data" (dark squares) combined with 300 kton-years of "antineutrino-beam data" (light 
[green] circles) for the following input value for the mixing parameters (indicated by the [blue] stars): Amfa = 3 x 10"'^ eV^, 



= 0.008, and sin^ 



= 0.5. See text for details. 



Several other comments can be made regarding Figs. |[ ^, and ^. First, as hinted in Sec. ||, the size of the sin^ 623 
uncertainty is significantly larger for maximal atmospheric mixing than for the nonmaximal cases. Second, even 
though sin^ 6*23 = 0.35 and sin^ ^23 = 0.65 yield approximately the same value of sin^ 26*2! and therefore similar 
number of disappearing v^s (see Sec. they yield very different number of appearing i/gS, as discussed in Sec. III. 
This fact is clearly visible if one compares Figs. |[ |[ The error contours are smaller for sin^ 623 = 0.65 (Fig. |^) given 
the larger statistics. 

In order to "solve" the degeneracies, more information is needed. One option is to run with an antineutrino beam 
and keep the same experimental setup, or perform a different experiment with a different baseline and neutrino beam. 
We will here pursue the first option, but warn readers that running with antineutrinos takes a significantly larger 
amount of running time in order to obtain the same amount of data (especially if the neutrino mass hierarchy is 
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FIG. 4: Same as Fig. H, for sin^ 023 = 0.35. 



normal). In our opinion, antincutrino running is only realistic if an improved proton driver, such as the one currently 
being investigated was present at NuMI. The existence of another experiment with a different neutrino beam and 
baseline might turn out to be realistic if the JHF to SuperKamiokande proposal materializes. 

We depict the two-sigma CL measurements in the (Amfg x |t/e3p), (|?7e3p x sin^ ^23), and (Ato^s x sin^ 023)-planes 
obtained after 120 kton-years of "neutrino-beam data" and 300 kton-years of "antineutrino-beam data" in Figs. ||, IJ, 
and ^. One can see that in the case of maximal atmospheric mixing (Fig. |) the two-fold de generacy is lifted, and 
not only does one obtain a better measurement of |J7e3p but, more importantly (in our opinion), one can determine 
the neutrino mass-hierarchy. In the case of non-maximal atmospheric mixing, the four-fold degeneracy collapses to a 
two-fold degeneracy, and the neutrino mass hierarchy can also be unambiguously measured. 

The capability of determining the neutrino mass hierarchy will depend on the size of |C/e3p. Fig. I depicts the 
value of Ax^ 
value of Am? 



IXmin 



X^ijj(Am^3 > 0)1 as a function of the input value of |?7e3p, assuming that the real 
Ax^ is defined after integrating out the marginal parameters |C/e3p and sin^6'23, 
and is to be interpreted as a Ax^ for one degree of freedom (Am^g).^ This means that, if sin^ 6*23 = 0.5, a two sigma 
hint that the neutrino mass hierarchy is normal should be obtained if |C/e3p = 0.003, while a five sigma discovery 
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(Am^g <0) 
3 X 10-3 eV^. 

2 



5 This is distinct from the analysis presented in [ 



8 



> 

CD 



E 

< 



0.005 
0.004 
0.003 
0.002 
0.001 


-0.001 
-0.002 
-0.003 
-0.004 
-0.005 



10 



0.3 0.4 



0.5 



0.6 0.7 



10 



-3 



_ 1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 _ 


_ 


1 1 


1 1 








1 


1 1 








_ 


- 


















(iitm 

m 














E 














P*llH* 








- 


























- 

- 


























: 








































































»iliHi> 


iHli 


I- 


~- 










1 1 
















" 1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 ~ 


1 1 








i 


i i 


















-3 -2 
10 10 ,,, ,2 


m 








mSSm 








~- - ~- - - 








■.•.■Kz:-.: 





























































1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 









0.005 
0.004 
0.003 
0.002 
0.001 


-0.001 
-0.002 
-0.003 
-0.004 
-0.005 



0.3 0.4 0.5 0.6 0.7 



sin^e 



23 



FIG. 5: Same as Fig. H, for sin^ 023 = 0.65. 



of the neutrino mass ordering should be achieved if |J7e3p ^ 0.009. 
of sin^ 023 , and increases with increasing sin^ 9 



The discriminatory power depends on the value 
23. As has been pointed out in Sec. Ill, this is due to the fact that 
the number of Ve events due to oscillations is proportional to sin^ 623 , and more statistics can be obtained after the 
same running time for larger values of sin^6'23. Note that this is true even though the solution with sin^ 6123 7^ 0.5 



(statistically) is still ambiguous when it comes to determining the value of |C/e3| whether 6*23 > 7r/4 or 623 < 7r/4. 



As mentioned in Sec. II], solving the 623 ^ tt/2 — 623 degeneracy is much trickier. It is not resolved by comparing 
the neutrino channel with the antineutrino channel, or even comparing different searches for —f i/g-transitions in 
long-baseline experiments, as long the matter effects are far from the a ~ 1 "resonance" region. In order to probe 
such a regime, however, one is forced to probe very large neutrino energies, and hence very long neutrino baselines, 
which implies the construction of new, futuristic, neutrino facilities. 

As we pointed out in Sec. Ill, other options include looking for i^7.-transitions [ p5| . In order to acomplish this, 
it is necessary to build an intense, high energy source of electron-type neutrinos, i.e., a neutrino factory. One would 
have to decide, however, if the construction of such a device is worthwhile given its physics capabilities in the scenario 



9 




FIG. 6: Ax^ = |Xmin(Amf3 < 0) - Xmin(Am?3 > 0)1 as a function of \Ue3f , for Amis = +0.003 eV^ and sin^ 6I23 = 
0.35, 0.5, 0.65, obtained after 120 kton-years of "neutrino beam data" and 300 kton-years of "antineutrino beam data." See text 
for details. The solid horizontal lines indicate the two- and three-sigma confidence levels for determining the sign of Amf^. 



we are considering here. " 

One final option would be to study the disappearance of electron-type (anti)neutrinos from nuclear reactors. Unlike, 
say, the KamLAND experiment, such reactor would not have to be a very long baseline setup, and CHOOZ-like 
distances (0(1 km)) would be optimal. The challenge would be, however, to improve the "sin^ 20-sensitivity" of the 
CHOOZ experiment by, say, a factor of 10, from 0(0. 1) to 0(0.01). Such an experiment would have to deal, for 
example, with the issues which help define the sensitivity of the CHOOZ experiment, such as the overall neutrino flux 
normalization, which is nominally known to a few percent |]l9| . One way to try and address this issue would be to 
also include a near detector, a la Bugey pGf . 

Next, we would like to address a more "practical" but important issue, related to how well should the atmospheric 
parameters be measured before a "safe" decision regarding the location of the off-axis detector can be made. Part of 
the information required in order to address this issue is already contained in Fig. One sees that the sensitivity 
worsens significantly for |Ato^| 2 x 10""^ eV^ because the oscillation maximum begins to "leave" the energy window 
we are looking into, and the statistics begins to deteriorate significantly. A different question we would like to address 
is whether one can still try to obtain information regarding the neutrino mass hierarchy. We do this by repeating the 
analysis depicted in Fig. || for fixed sin^ 6*23 = 0.5 but significantly different values of Ato^3. Fig. |^ depicts the value 
of Ax^ as a function of the input value of |?7e3p for Amfg = -|-(2, 3,4) x 10"'^ eV^, which are meant to cover the 
currently allowed 99% CL measurement of jAm^gj. A more precise measurement should be expected after a few more 
years of K2K data ||l|, and will definitely be obtained after MINOS (22| and the CNGS[|3| experiments "turn on." 

One should note that, even for Am^g = 0.002 eV^, where the statistics is expected to be significantly worse, a three 
sigma effect can be obtained for |/7e3p ^ 0.01. For Amjg = 0.004 eV^ it turns out that the discriminatory power is 
larger than for Amfg = 0.003 eV^, even though the statistical sample is slightly smaller. This is due to the fact that 
for Amfg = 0.004 eV^ the maximum of the oscillation takes place at a larger value of the neutrino energy (compared 
with the "optimal" Amfg = 0.003 eV^ ), such that when the mass hierarchy is inverted, the oscillation maximum 
starts to "leave" the energy window we are looking into. This means that the number of events obtained with the other 
choice neutrino mass-hierarchy is "more different" than in the case Am^g — 0.003 eV^, improving the discriminatory 
power. Incidently, the same phenomenon renders the discriminatory power in the case Am^g = 0.002 eV^ worse. Note 
that this effect came about because we chose the input value of Am^g to be positive. Were Amfg negative, we would 



^ We do not advocate that there is no case for building a neutrino factory if the LMA solution is ruled out. We do believe, however, that 
its physics goals would have to be reviewed in detail. This statement, of course, only applies if no other new physics is discovered in the 
leptonic sector untill then. 
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have observed the opposite trend. 

Finally, one may also inquire what happens if the baseline were different from the 900 km we have chosen here. It 
turns out that for distances which are slightly smaller than 900 km {L > 700 km) similar results can be obtained, 
as was discussed in detail in [|l^. For significantly smaller baselines {L ^ 500 km), however, the matter effects 
become less significant, and the "optimal" energy for a fixed value of Amf^ becomes significantly smaller. This 
implies that while one might still keep the sensitivity to observing |?7e3p to similar levels, the power to discriminate 
between the two neutrino mass-hierarchies disappears. Furthermore, it should be added that performance for the 
iron-scintillator detector we are assuming for neutrino energies between 1 and 3 GeV does not apply for an equivalent 
energy window centered below 1 GeV, as the capability to distinguish electrons from neutral current events (and even 
muons) deteriorates. In this sub-GeV regime, one might be better off using different detector technologies, such as 
Water Cherenkov detectors (as is being proposed for the JHK to SuperKamiokande project)** 



V. CONCLUDING REMARKS 



In the most conservative of scenarios, the fact that neutrinos have mass and mix will be confirmed as the solution 
to the atmospheric and solar neutrino puzzles, while the LSND result will be ruled out as a signal for neutrino 
oscillations. In such a world, the main goal of the next generation of neutrino oscillation experiments is very clear: 
measure (or further constrain) the "small" Ue3 element of the leptonic mixing matrix. This is true irrespective of 
what the solution to the solar neutrino puzzle is. If the solar mass-squared difference is large enough, there is also 
the exciting "bonus" of trying to probe CP-violation in the leptonic sector. 

If the solution to the solar neutrino puzzle requires a very small Am^^j^j., CP violation becomes inaccessible, as 
"solar effects" simply cannot be observed in terrestrial experiments. Nonetheless, we have shown that the physics 
prospects for future long-baseline neutrino experiments are still very attractive. In many respects, measuring Ue3 
becomes a simpler task, while several obstacles still remain. While overcoming these obstacles, we have shown that 
one can (must) not only measure |C/e3p in long-baseline experiments with non-negligible matter effects, but also 
determine the neutrino mass-hierarchy. Determining the mass hierarchy is arguably the second clear goal of future 
neutrino experiments (whether or not neutrino oscillations are directly involved) and just as exciting (in our opinion), 
as discovering CP-violation in the leptonic sector. Whether or not the neutrino mass-hierarchy is normal or inverted 
is certainly cleaner to address if "solar effects" are not around to complicate things, and requires information from two 



* We thank Michal Szleper for comments concerning detector performance issues. 
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distinct oscillation channels, such as — s- i^e-transitions at a long baseline, complemented by — > j/g-transitions at 
the same experiment. Another option is to combine a "long" and a "short" long-baseline search for — > i^e-transitions. 
This might well be the scenario that we will contemplate in the future if both the JHF to SuperKamiokande and the 
NuMI off-axis efforts materialize. 

We have included the effects of "translating" the information obtained in the disappearance channel to the appear- 
ance channel, error bars included. Not only can we cleanly identify the confusion induced by not knowing what the 
neutrino mass-hierarchy is, but also identify the ^23 versus 7r/2 — ^23 degeneracy [ p^ , which proves to be much harder 
to lift in the type of long-baseline experiments we (and most of the community) imagine as feasible in the near to 
intermediate future. We do provide a few candidate solutions, which involve new neutrino oscillation experiments, 
such as a "Super-Bugey" effort |^ or hunting for transitions. Finally, it should be noted, of course, that 

several of the issues we raise here are also applicable if the LMA solution is confirmed by KamLAND. 

In summary, the results of the KamLAND experiment (and, in a few years, of MiniBoone ||2^) will shape the 
future of neutrino physics experiments. If we are on the right track (i.e., if there are no major surprises or set- 
backs), however, all different roads seem to lead to a next-generation, long-baseline experiment capable of searching 
for Vp, i^e -oscillations. Even if the LMA solution is ruled out, several fundamental measurements can be performed 
by such experiments, which are worthy of very serious consideration. 
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APPENDIX A: OSCILLATION PROBABILITIES AND THE LOW SOLUTION 

The effective Hamiltonian that describes the time evolution of neutrinos in matter can be written in the flavor basis 



as 

2E 



^13 







J c/t + 









H=^{U\ Ami, \U^+\ ]}, (A.l) 



where a — G fV^N^^E represents matter effects from the effective potential of electron-neutrinos with electrons, and 
U is the flavor mixing matrix in vacuum (PDG representation), 



C12C13 S12C13 sise 



U = \ -S12C23 - Ci2S23Sl3e** C12C23 - Si2S23Sl3e*'^ S23C13 | . (A. 2) 

S12S23 - 2120235136"^ -C12S23 - Sl2C23Sl3e*'' C23C13 

By imposing jArTi^gl > Amfj > 0, we can relate ^23 and jAm^gl with the atmospheric angle and mass-squared 
difference, 6*12 and Ami, with the solar angle and mass-square difference, and ^13 with the "reactor" angle. These 
definitions are used throughout the paper (see (l6j for details). 
For a baseline L, the evolution of neutrino states is given by 

i^{L) = S{L)iy{0), (A.3) 

with 

S{L) = e""^ (A.4) 

for constant matter density. The corresponding effective Hamiltonian for antineutrinos is obtained by U U* and 
a —a. 

If the Hamiltonian can be separated into H = Hq + eHi, where Hq can be exactly solved and e/lll such that Hi 
can be appropriately treated in perturbation theory, then S{L) — Sq{L) + Si{L) and So{L) = e~'^oi gives the lowest 
order transition amplitudes. 
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For lAm^jl ^ Irfm^al, and non negligible a, we choose 











Ami 



Ami 



(A.5) 



where Am^ describes the energy-level spacing in matter and U the effective mixing elements, and ignore Hi effects, 
which are suppressed by e = Amf2/Amf3. Of course, 



So{L) = U 



V 



(A.6) 



The exact diagonalization of Hq yields the effective mass differences and mixing angles in matter in the limit 
\Ami2/ Amf^l 0, while no approximation have to be made regarding a. One obtains 



U : 



( -—{h 

-C23 



-13; 



-13 



-^Si-iSx-iCx-i — S23S13C13 
;77C23Si3Ci3 ^ 



\ S23 ;7;-C23Sl3Ci3 — C23S13C13 / 

for the mixing angles in terms of the vacuum parameters, where we have defined 



and 



'3 = 



,2 

13' 



(A.7) 



(Ai 



Ato?3 
Am% 
Amf3 



1 + a - -v/l + a2 - 2acos(26li3) 



1 + a + a/I + - 2acos(26'i3) 



(A.9) 



for the sake of 



which give the effective mass differences. We have introduced the dimensionless parameter a = 
simplicity. The ordering of levels in matter is (1,2,3) for the hierarchical case of Am^3 > and a < 1. 

The effective mixing matrix U is independent of 612 and 6, if the latter is appropriately rotated away. The matrix 
U can be expressed also in the PDG form. To do so, we change 



1 



(A.10) 



By comparing with previous expressions, we may establish the correspondence between the effective mixing angles, 
6ij, and the vacuum parameters, getting 



C12 = 0, 



Cl3 



C23 = C23, 
I2 



-13 



n2 



Sl3 = 



|5l2| = 1, 

S23 = S23, 

2 
13 



"3 



(A.11) 



up to signs. The vanishing mixing in matter C12 = is a consequence of the degeneracy Am^2 = in vacuum and 
says that the lowest mass eigenstate in matter contains no electron- neutrino flavor component. This result is the 
ingredient that avoids genuine CP violation in matter, even if one has three non-degenerate effective masses. 

The transition amplitudes for I'/s in the Am^2/Am^l3 — > case are given by matrix elements, which can 

be written 



A{a^f3;L) = 5o(L)/3a 



(A.12) 
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From this expression we may calculate all probabilities, 



where 



and 



P{^^. 



■ Ve) 



1 - sin^(2^i3)sin^ 
sin^ (26113) sin^ 



Al3 
Ai3 
Ai3 



C23 sin^(2^i3) sin^ 
1 - 5^3 sin2(2^i3)sin2 
cos(2^i3)sin [Ai3(l-f 



Ai3 

a)] sin 



2S23C23 {1 ~ [Ai3(l + a)] cos 



Ai 



Ai 



} 



2 2 
'*23'^23 



|2-2cos[Ai3(l + a)]cos A13 - sin^ (2^13) sin 



A 



13 



2 cos(26ii3) sin [Ai3(l + a)] sin 



Ai3 
2 „2 



} 



P{vr^yT) = 1 - C23sin^(26i3)sin^ A13 - 2S23C23 |l - cos [Ai3(l + a)] cos 



A 



13 



cos(26'i3) sin [Ai3(l + a)] sin 



^13 



} 



^13 



Ai3vr+ 



i;2-2acos(26'i3), with A13 = . 



sin2(20i3) =4 



„2 2 
*13''13 



1 + ~ 2acos(26'i3)' 



(A.13) 
(A. 14) 
(A.15) 

(A.16) 

(A.17) 

(A.18) 

(A.19) 

(A.20) 
a. The 



In order to get the corresponding expressions for antineutrinos, one simply exchanges a — > —a, i.e., a - 
effect of such a change in the probabilities comes from the different relative sign between mass and matter terms in 
Hq. It is important to note that, because we are disregarding Ami2 effects (which is perfectly justified, as we are 
interested in the LOW solution) the expressions for neutrinos and antineutrinos are exchanged when the neutrino 
mass hierarchy is "flipped", i.e.. 



(A.21) 



This behavior is a consequence of the fact that we are, in practice, dealing with an effective two-level system. 
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